function out=vech(in)
nxx=size(in,1) ;
Lc=v2h(nxx) ;
out=Lc*vec(in) ;
end
